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Radar Backscattering by The Icy Galilean Satellites and The Radar Glory Effect. 
Eric M. Gurrola and Von R. Eshleman. Center for Radar Astronomy, Stanford University. 


Radar echoes from the icy Galilean satellites — Europa, Ganymede, and Callisto — are unique 
in character, and they may be explained by a phenomenon that has been called the radar glory 
effect. 1 ’ 2,3,4 We have analyzed the radar glory backscattering from buried craters as a possible 
model for this effect. These craters have a smaller refractive index below the crater than above 
as shown in Fig. 1(A). The possibility exists that the rays shown will be totally internally 
reflected at the crater walls, and we assume this to be true in the following. The rays which 
contribute to the backscattering will come from a circular annulus, when viewed before the 
refraction occurs at the ice-vacuum surface, in a plane orthogonal to the rays reflected from 
the crater (see Fig. 1(B)). For craters near the limb of the satellite some of the rays may be 
blocked by the crater rim, and this circular annulus will break up into two circular arcs. We 
assume that these glory arcs are broken up into glints of dimension H , assumed larger than the 
wavelength, and that the fields from a pair of glints connected by the ray shown add coherently 
while the fields from the various glint pairs add incoherently. The coherent addition of the fields 
from glint pairs arises because for the ray shown in Fig. 1(A), there is another ray that travels 
the exact same path, but in the opposite direction. 



Figure 1. (A) Geometry of a buried crater and a three-bounce (N=3) ray. The angle 
0 is between the surface normal of the satellite and the radar line of sight. For an 
ice-overburden, n 0J = 1.8. In the model n 0 2 < n 0 i (i.e., n<l). (B) Glory arcs of angular 
extent 2<fo and a pair of glints of dimension H . TE and TM electric field vectors are 
also shown. 


Consider a crater located near the surface of the satellite at coordinates (/2, 0, <£), where R 
is the satellite radius, 0 is longitude with 0 = 0 in the direction of Earth, and $ is latitude with 
$ = tt/ 2 in the plane defined by the spin axis of the satellite and the 0 = 0 direction. The radar 
cross-sections for fields received in the orthogonal and same linear (OL and SL respectively) 
polarizations and those for fields received in the orthogonal and same circular (OC and SC 
respectively) polarizations are 


^sl(Q,3>) = CA L (3 - X L ) , CT OL (0,$) = CA L ( 1 + X h ) (1) 

0sc(Q) = 2CA c (l + X c ) , a oc (Q) = 2CA c (l-X c ). (2) 
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The constant C is proportional to the number of glint pairs in the glory arcs, and is given by 
C = irrlFH/\ where F is the fraction of the glory arcs filled by glints, H is defined in Fig. 
1(B), and A is the radar wavelength. The A ' s and X’s are given by 

A h = + Tl) + (2? - T 2 2 ) cos2($ - i P )][<h - (1 + *)«] + + x)S) (3) 


A l X l = do{ 2 TiT 2 x</>o 

+[(T X 2 - Tf) cos2(* - i P ) - (T 2 + T 2 2 )cos4(* - */»)][*> - (1 + *)«] (4) 

+27iT 2 [3(1 + x)6 - x<fo]cos4($ - ip)} 

Ac = d 0 {(r 2 + T 2 )<fo - (T x - T 2 ) 2 (l + x)6} (5) 

A C X C = 2 aoTiTixfo (6) 


where 


. _ sin 26 cos 2 0 
a ° N nj) X cos 2 0' ’ 


(7) 


tiqi sin 0' = sin0, ip is the direction of polarization for the linear polarization case, T\ — 
4noi cos 0 cos ©'/(cos 0 + «oi cos 0') 2 is the product of the TE Fresnel transmission coefficient 
for the ray entering the ice and the TE Fresnel transmission coefficient the ray exiting the 
ice, T 2 = 4n 0 i cos 0 cos 0'/(noi cos 0 + cos0') 2 is the same product for the TM fields, x = 
cos [2 TV tan -1 (cos 0(sin 2 0 - n 2 ) 1 / 2 / sin 2 (9)] is the mode decoupling factor defined in reference 1, 
6 = </>o/4 — sin(4<fo)/16, and = w/2 if no shadowing occurs or sin <p o = cot 0 cot 0' when this 
last definition has a real solution. 

The power spectra are sums of the contributions from several craters along lines of constant 
Doppler shift which are lines, when projected onto the plane of the sky, orthogonal to the line 
of sight and parallel to the spin axis. That is, the cross-sections at a specified frequency are 
integrals of the cross-sections in equations (1) and (2) on semi-circles on the satellite surface 
where cos 0 cos $ =constant. In Fig. 2 we compare the computed cross-section for circular 
polarization with echoes obtained from Eurpoa; this spectrum is typical of the spectra obtained. 

The total cross-section of each satellite is an integration of the cross-sections in equations (1) 
and (2) over the entire illuminated surface. These total cross-sections can be written in a form 
similar to those in equations (1) and (2) except that the A’s and X’s are the same for linear 
and circular polarizations and are complicated integrals that must be computed numerically. 
In Table 1 we compare a summary of the main features of the data with the model predictions. 
The scaling constant C for these total cross-sections is given by C = {M irr^ R 2 )(F I! / X) 
where M is the total number of craters on the surface. 
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Figure 2. OC, SC, and OC+SC cross-sections obtained from Europa (dashed curves) 
on January 10, 1990 at the Arecibo Observatory at 12.6 cm-A and from the model (solid 
curves) using N — 3, n=0.7967. The data and the model have been scaled to give a 
total cross-section of 2.6. 



Europa 

Ganymede 

Callisto 

Buried craters 
N = 3, n=0.7967, C = 1 

b 

2.60 

1.52 

0.64 

0.67 

bsc/ doc 

1.56 

1.55 

1.19 

1.60 

^ol/^sl 

0.47 

0.47 

0.55 

0.44 

m 

1.73 

1.46 

1.43 

1.95 


Table 1: The data for the Galilean satellites are for 12.6 cm-A and are from Ostro 5 . 
b = T d 0 c = d$L + d 0 L- These cross-sections are normalized by the geometrical 
area irR 2 , i.e., b = a / tcR ? . Here m is the exponent for a least squares fit of the power 
spectra to a cos m 0 law. 
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